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\ Abstract. Studying the order of the chkal transition for Nf =2 is of fundamental importance to 

understand the mechanism of color confinement. We present results of a numerical investigation on 
the order of the transition by use of a novel strategy in finite size scaling analysis. The specific heat 
) , and a number of susceptibilities are compared with the possible critical behaviours. A second order 

^— ( ■ transition in the 0(4) and (9(2) universality classes are excluded. Substantial evidence emerges for 

' a first order transition. Results are in agreement with those found by studying the scaling properties 

of a disorder parameter related to the dual superconductivity mechanism of color confinement. 
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INTRODUCTION 



^ . Experiments indicate that confinement is an absolute property of matter. Indeed the 

upper limit on the number of free quarks per proton is R = riq/np < 10^^^, while 
rin ! 10^^^ is expected from the Standard Cosmological Model. A reduction factor 10^^^ 

p • is difficult to explain in natural ways unless R = 0, which means that confinement is 

1^ . related to some symmetry of the QCD vacuum. This also implies that the deconfining 

transition is associated to a change of symmetry of the vacuum, i.e. it is an order-disorder 
transition. This scenario must be confronted with direct studies of the deconfining 
^ ■ phase transition: experimental studies going on with heavy ion collisions have not given 

definite answers yet and most of our knowledge relies on numerical simulations of QCD 
on the lattice. From this point of view, the case of two light degenerate dynamical 
flavours is of special interest. A schematic view of the phase diagram for Nf = 2 is 
shown in Fig.[T] m is the quark mass and jU is the baryon chemical potential. 

In the /i = plane, as m oo quarks decouple and the system tends to the quenched 
limit. There the deconfining transition is an order-disorder first order phase transition, 
Z3 is an associated symmetry and the Polyakov line (L) is an order parameter. Z3 is 
explicitely broken by the inclusion of dynamical quarks and (L) is not a good order 
parameter, even if it works as such for quarks masses down to m ~ 2.5 — 3 GeV. 

At m ~ there is chiral phase transition at Tc ~ 170 MeV, from a low temperature 
phase where chiral symmetry is spontaneously broken to an high temperature phase in 
which it is restored: the corresponding order parameter is the chiral condensate {xf^Xj^) . At 
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FIGURE 1. Schematic phase diagram of Nf = 2 QCD. 



some temperature Ta > Tc also the Ua{^) symmetry, broken by the anomaly, is expected 
to be effectively restored. It is not clear which relation exists between the chiral transition 
and the deconfining transition: empirically the Polyakov line has a rapid increase at the 
transition temperature, indicating deconfinement. The transition line in Fig. [T] is defined 
by the maxima of a number of susceptibilities {Cy, Xm, ■ ■ all coinciding within errors, 
which indicate a rapid variation of the corresponding parameters across the line. 

At m ~ a renormalization group analysis plus e-expansion techniques can be made, 
assuming that the relevant degrees of freedom for the chiral transition are scalar and 

^ni 

pseudoscalar fields [1, 2, 3]. If the Ua{1) symmetry is effectively restored, i.e. if the rj 
mass vanishes at Tc, then there is no IR stable fixed point and the phase transition is first 
order; if not an IR fixed point exists, which can produce a second order phase transition 
in the 0(4) universality class. 

In the first case the transition is first order also at m 7^ and most likely up to m = 00. In 
the second case a phase transition is only present at m = 0, which goes into a continuous 
crossover as m 7^ 0: that means that one can move continuously from confined to 
deconfined and that no true order parameter exists. This would be in contradiction with 
the deconfinement transition being associated to a change of symmetry: the issue is 
therefore fundamental. 

The problem has been investigated on the lattice by several groups with staggered |0, 
or Wilson [11] fermions. The strategy used has either been to search 
for signs of discontinuity at the transition, or to study the scaling with m of different 
susceptibilities, or to study the magnetic equation of state. No clear discontinuities have 
been observed, but also no conclusive agreement of scaling with 0(4) critical indexes. 
We present the results of a big numerical effort aimed at clarifying the issue. We use 
non improved Kogut-Susskind action and lattices 4 x with = 16, 20,24, 32. A full 
account of our results can be found in Iil2i1 . 
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FIGURE 2. Specific heat peak value for Runl (left) and for Run2 (right), divided by the appropriate 
powers of to give a constant. Both the 0(4) and 0(2) critical behaviors are displayed. 



RESULTS 

The theoretical tool to investigate the order of a phase transition is finite size scaling. 
The extrapolation from finite size to the thermodynamical limit is governed by critical 
indexes, which identify the order and the universality class of the transition. Around the 
chiral transition the system has two fundamental lengths: the correlation length ^ and the 
inverse quark mass l/mq. ^ is usually traded with the reduced temperature T = l—T/Tc, 
^ ~ T^^ as T ^ 0. The effective action depends on the order parameter, as dictated by the 
symmetry, and as t ^ irrelevant terms can be neglected; the correlators of the order 
parameter describe the thermodynamics. The most important quantity is the specific 
heat, which shows the correct critical behaviour independently of the identification of 
the order parameter. 

For the specific heat and the susceptibility of the order parameter the scaling laws are 

Cv - Co ~ Lf />, (^zLy^am^Lf^ ; (1) 
Z^LJ/>;^(TLy^am,L^'■) . (2) 

Co stems from an additive renormalization. An alternative way to write them is 

Cv-Co- L'^'^c (T(am^)-i/(^-^''),am^L^'') (3) 
X^Lj/'h {T{am,)-'/^''y'^\am,Ly/^ . (4) 

The values of the indexes characterize the transition: the values relevant to our analysis 
are listed in TableHl 0(4) is the symmetry expected if the transition is second order, but 
it can break down to 0(2) by lattice discretization for Kogut-Susskind fermions ITJ at 
non zero lattice spacing. 

The scaling analysis is made difficult by the presence of two independent scales. The 
attitude taken in the previous literature has been to assume the volume large enough so to 
neglect the dependence on Li. since at fixed amq, j5 the susceptibilities must be analytic 
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FIGURE 3. The same as figure|2]for the chiral susceptibiUty Xm- 



in the thermodynamical limit, at large Ls the dependence on aniqlZ^ must cancel the 
dependence on Ls in front of the scaling functions in Eq.s ([T]) and ©. It follows that 
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The peaks of (Cy — Cq) and of x should then scale as 
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as aniq —>■ 0. The positions of the maxima scale as x{amq) V(v>7,) = const. 

One can also consider to keep xL^J^ fixed while taking aLs ^ l/m^. This assumption 
should work better if is still comparable to the correlation length, which may be the 
case close enough to the critical point. In this case the scaling laws are 
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const and the 



Eq.s dVl) stay unchanged, the positions of the maxima scale as tL, 
width of the peaks are volume dependent. 

We have instead followed a novel strategy which does not rely on any assumption: we 
have performed a number of simulations at different values of Ls and aniq keeping the 
variable aniqlJs'' fixed and studying the scaling with the respect to the other one. In doing 



TABLE 1. Critical exponents. 
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so one has to assume a value for y^: we have chosen that expected for 0(4), which is the 
same within errors as for 0(2). From Eq.s dU) and Q it follows that, at fixed aniqLl'' 

{Cv - Co)max °^ Lf^^ ; Xmax °^ LJ^^ ■ (9) 

as ^ oo. If 0(4) or 0(2) is the correct symmetry, the values of a/v and 7/ v should be 
consistent with the corresponding values listed in Table[T] We have run two such sets of 
Monte Carlo simulations, called in the following Runl and Run2, with aniqlJj' = 74.7 
and aniqlZ'' = 149.4 respectively. The spatial lattice sizes used for each of the two 
sets are L5 = 12,16, 20, 32, the standard hybrid R algorithm 11 131] has been used to update 
configurations. In Figs.|2land|3]we show the peak values of the specific heat and of the 
chiral susceptibility, divided by the appropriate power of L^ , as a function of Ls (see 
Eq.|51): scaling is clearly violated, 0(4) and 0(2) universality classes are excluded. 
An alternative way to study the order of a transition is to look at scaling of pseu- 

docritical couplings: one can try the two alternative scaling laws Tc = k^L^J^ or Xc = 
fc!^(am^)^/(^^'''^ . The physical temperature T = \ / {Lta{^ ,mq)) is a function of both /3 
and aniq, so that the reduced temperature T can be expanded as a power series in (/3 — /3o) 
and in aniq, where /3o is the chiral critical coupling. Only the linear term in /3 was con- 
sidered in the previous literature. We have found that the following terms are sufficient 
to fit the data 

'^'^ iPo -P)+kinamq + k,„2{amq)^ + k^pamq{l5o-P). (10) 

Our result is that it is not possible, within the present mass range, to discriminate 
among the possible critical behaviours by looking at pseudocritical couplings only. The 
inclusion of other terms in Eq. (flOt . besides the one linear in /3 solely considered in 
previous literature, is however crucial to obtain any scaling at all. 

We can also use our data to perform the scaling analysis in the same way as done 
in previous literature, i.e. by making some assumption on the reaching of the thermo- 
dynamical limit: no universality class is chosen a priori in that case and one can test 
all the possible critical behaviours. We have found that assuming that aL^ 3> l/mji but 
still Ls ^ ^, i.e. Eqs. dU), works better: this is reasonable when ^ goes large. We have 
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FIGURE 5. Scaling of the equation of state for (9(4) (left) and first order (right). 



added to the data from Runl and Run2, those from two other simulations performed 
at Ls = 16, aniq = 0.01335 and Ls = 24, aniq = 0.04444. In FigHwe show the scaling 
obtained for the specific heat: 0(4) is again clearly excluded, while a good agreement is 
found with a weak first order critical behaviour. A similar behaviour is observed for the 
chiral susceptibility. 

We have also studied the magnetic equation of state, i.e. the scaling of (v^V^) — 

{V/V/)o = aml^^Fizml^"^''). Results are shown in Fig. 13 again the first order behaviour 
describes well the data while the second order is excluded. 

Even if we have still not found clear signs of discontinuities in physical observables, 
our results give some evidence for a weak first order chiral transition which would 
persist also at m 7^ 0. That would be in agreement with the deconfining transition being 
order-disorder and with confinement being associated to some symmetry of the QCD 
vacuum. This naturally leads to look for such a symmetry and for an order parameter 
which describes the deconfining transition both with and without dynamical quarks. 
One candidate is dual superconductivity of the vacuum and the related parameter is the 
vacuum expectation value of a magnetically charged operator (/i). Indeed it has been 
shown that it is good order parameter also in full QCD [14] and its critical behaviour 
across the Nf = 2 transition has been studied in Ref. 1 15]. In Fig.|^we show the scaling 
of its susceptibility p = ln(/i): again results seem to indicate a first order critical 
behaviour. 



CONCLUSIONS 

We have argued that the study of the order of the chiral phase transition for Nf = 2 is of 
fundamental importance to understand confinement. We have shown that the analysis of 
pseudocritical coupling alone cannot discern between the possible critical behaviours. 
By adopting a novel strategy which reduces the finite size scaling analysis to a one scale 
problem, we have been able to exclude a 0(4) (0(2)) second order critical behaviour, 
while we have found consistency with a weak first order critical behaviour both in the 
scaling of susceptibilities and in the equation of state. This would be in agreement with 
confinement being an absolute property of matter related to some symmetry and with the 
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FIGURE 6. Scaling of the susceptibility p of the disorder parameter (jj.) for (9(4) (left) and first order 
(right). 



deconfining transition being order-disorder: indeed the analysis of the scaling properties 
of the susceptibility of an order parameter associated with the dual superconductivity 
mechanism of color confinement leads to analogous results. However we have still not 
found any clear evidence for discontinuities in physical observables. The issue is still 
open and we plan in the future to investigate it more deeply by making simulations with 
improved actions and algorithms and with aniqLl'' fixed according to first order. 
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